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THOMSON  BACKSCATTERED  X-RAYS  FROM  AN  INTENSE  LASER  BEAM 

I.  Introduction 

Tunable,  near  monochromatic,  high  brightness  x-rays  would  be  an  important  tool  in 
research  and  medical  diagnostics.  Synchrotron  light  sources  have  produced  useful  x-rays 
for  a  large  user  community.  In  this  paper  we  examine  a  closely  related  method  of  x-ray 
generation,  Thomson  backscattering  of  x-rays  from  intense  laser  beams. The  schematic 
is  shown  in  Fig.  1,  where  an  electron  beam  intersects  an  incoming  laser  pulse.  Radiation  is 
backscattered  at  a  double  Doppler  upshifted  frequency.  The  laser  pulse  in  Fig.  1  acts  in  a 
similar  fashion  as  the  static  magnetic  wiggler  in  synchrotron  light  sources  or  free  electron 
lasers.^~^^ 

One  advantage  to  this  approach  is  that,  because  the  wavelength  of  the  laser  is  many 
orders  of  magnitude  smaller  than  that  of  static  undulators,  an  electron  beam  of  much  lower 
energy  can  be  used  to  generate  x-rays  of  a  particular  energy.  For  example,  radiation  of 
0.04  nm  wavelength  (30  keV)  x-rays  can  be  generated  by  a  laser  with  1  fim  wavelength 
and  electron  beam  energy  of  40  MeV. 

An  experiment  is  under  way  at  Vanderbilt  University  to  utilize  the  FEL  as  the  undn- 
lator.  The  goal  is  to  develop  diagnostic  medical  imaging  techniques  based  on  the  detection 
of  atoQUc  species  importsuat  in  biological  substimces  by  utilizing  the  discrete  K-edges  which 
are  in  the  low  keV  energy  range.^  A  Isiser  undulator  has  also  been  proposed  as  the  damping 
mechanism  of  a  very  low  energy  (~  1  MeV)  storage  ring^*  and  for  emittance  reduction.** 
The  analytical  expression  of  the  Thomson  backscattered  radiation  from  a  laser  undulator 
derived  in  this  paper  can  be  applied  to  all  these  applications. 

II.  Formulation 

We  are  interested  in  calculating  the  Thomson  radiation  intensity  pattern  and  spec¬ 
trum  which  results  when  a  laser  pulse  intersects  with  an  electron  beam  with  small  initial 
transverse  momentum  going  in  the  opposite  direction. 

The  laser  pulse  is  assumed  to  be  linearly  polarized  with  frequency  ufi,.  The  vector 
potential  of  the  laser  pulse  can  be  separated  into  fast  and  slow  components, 

A(q)  =  A(»7)sin(T/)e*,  (1) 
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where  rj  =  ki,z  +  <*;£,/,  The  pulse  shape  A{rj)  and  the  wavenumber  ki,{Ti) 

are  assumed  to  be  a  constant  for  the  interaction  time  T  and  sin(77)  is  a  fast  oscillating 
component. 

The  energy  radiated  per  unit  solid  angle  (dfl)  per  unit  frequency  {du>)  per  electron  is 

|2 

,  (2) 


cPi 


rT/2 


dwdU 


/di  h  X  (n  X  ^)exp  [ia>(<  —  h  •  f/c)] 
-TI2 


where  n  is  a  unit  vector  pointing  from  the  radiating  electron  to  the  observation  point,  r 
is  the  electron’s  coordinate,  T  =  L/c  is  the  interaction  time,  L  is  the  length  of  the  laser 
pulse  and  (3  =  c~^dT/di  is  the  electron  velocity  normalized  to  the  speed  of  light  c.  We 
consider  an  electron  with  small  initial  transverse  velocity,  fixo  =  Wso/c  and  /3yo  =  Vyo/c. 
The  symbol  above  the  variables  denotes  a  function  of  (<). 

It  is  convenient  to  use  Cartesian  coordinates  for  the  velocity  and  position  of  the 
electrons  and  spherical  coordinates  for  the  Thomson  backscattering  radiation.* 


and 


n  X  (n  X  ^)  =  —  cos  6  cos  <l>  +  0y  cos  6  sin  <f>  —  sin  6)eo 
+  (4*  sin  4>  —  0y  cos 

n  •  f  =  X  sin  0  cos  ^  +  y  sin  ^  sin  ^  +  z  cos  B, 


(3) 


(4) 


where  =  ^  •  e*,  x  =  f  •  e*,  €uid  similarly  for  the  y  and  z-components. 
We  separate  the  two  components  of  the  radiation 


dwdU  dujdU  dufdn  ’ 

and  perform  the  calculation  in  the  variable  y,  where 


(5) 


dPL 


0 


dufdCl 


47r*c* 


/  dv 

J  —  Ari/2 


dx 


dy 


dz 


dr]  —  cos  B  cos  ^  cos  5  sin  ^  —  —  sin  B 


drj 


eW 


dufdSl  dTT^C* 


J-Av/2 


drj 


df]  or) 


.  ,  dy  ,\  r  T, 

Sin  0  -  —  COS  I  exp  1*V’] 


exp  [t^] 


(6) 


(7) 


cfidt  =  {dT/dT})dT],  At;  =  a>£,r  =  2irNo  and  No  is  the  number  of  periods  in  the  laser  pulse. 
In  the  1-D  limit,  there  are  two  constants  of  motion: 


Px  -  -A  =  P,o, 

c 


7(1-/S.)  =  7o(1-/9*o),  (10) 

where  P*  is  the  canonical  momentum  of  the  electron  in  the  x-direction,  and  the  subsci^ipt 
0  denotes  the  initial  value.^*~^* 

In  the  following,  we  assume  that  the  electron  transverse  motion  is  small,  i.e.,  a(^|o  + 
0lo)'yl/^  «  It  intensity  not  exceedingly  large,  i.e.,  a  «  27o//3*o»  and  fcj.Ar  <<  1, 
where  Ar  is  the  radius  of  the  electron  oscillation  in  the  transverse  direction  driven  by  the 
laser,  a  =  (e/miiC^)A  and  mo  is  the  rest  mass  of  the  electron.  The  particle  motions  are 


dij  {uil/c) 


0to - sinr;  , 

70  J 


(11a) 


dr]  {u>l/c) 


dz  1 
—  - 

drj  {ul/c) 


(fi)'- 


1  +  cos(2t7))  +  /?,o  —  8in(T;)  1  , 


(116) 


(11c) 


where  =  (1  -  7o  *  -  initial  axial  velocity,  =  (1  -  (1  +  7^(/3|o  + 

^yo))/(7o(l  +^*o)^))/2,  0x0  =  0zo/{^  +  0zo)i  0yo  =  ^yo/(l  + /?*o )  and  a  =  a/(l  +/3*o)- 


Their  locations  are 


(12a) 


y  =  7 - {v  +  At;/2), 

(wt/c) 


(126) 


(12c) 
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The  9  and  (f)  components  of  the  radiation  are  now 


d^Ie 

dudQ 


c2u,2  [■/_  _  \  \ 

^  -r-z — 5-  /  dr]  {  /9*o  cos  6  cos  <j)  +  cos  9s\fx<f>—  [  -  — -r  )  sin  5  ) 

47r2ca;*  |y_A,/2  [V  \  5  /  J 

5  -  .  .  1  ■  ^ 

- (cosScos^  + /3*o  sin^)sinT; - 2sin^cos(277)  exp[t^]  , 


d^I  f^v/2  a  1  - 

T  — 2  /  {0tosm<f>~  of^os<i>) - sin  sin  7/  exp  [itA]  ,  (14) 

awdil  47r  cufj^  J-Ati/2  To  J 


where 


=  V’o  H - [1  -  ()9*o6in5cos^  +  /Sj,o8intfsin^  +  ;9i(l  +  cos 5))]  {rj  +  ^rj/2) 

- (sin  6  cos  ^  +  /Sjco  ( 1  +  cos  ^))  /  drj'  —  sin  77' 

J—A71/2  To 

+  — (1  +  cosd)  f  7~2^v‘ 

J-Af,/2  4To 

07  n  .  . 

- (l  +  COs5)  /  — =•  COs(277  )d77' 

•^L  J-Av/2^1fo 


is  the  phase  and  ijjo  =  — (o7/a>£,)A7;/2  —  (w/c)zo(l  +  cos^). 

The  integral  in  Eq.  (15)  may  be  evaluated  to  obtain 

^  ~  V»o  +  dftij  +  dj.  cos  77  +  d,  5in(277), 


where 


(/9*o  cos  <j>  +  0yo  sin  4>)  sin  9 


-('-5) 


(1  +  cos^) 


U7  r  —  1  ** 

d*  =  —  sin<?cos^  +  /3*o(l  +  cos^)  — , 

To 


(17a) 


(176) 


d,  = - (1  +  cos9)^. 

07L  87$ 


Substituting  (16)  into  (13)  and  (14),  we  obtain 


d^Iff  ^  c^oj^ 

dci;dn  ~ 


2  9o,0lo - h  {coB9cos<f>  +  0ro  si"^)  -  7-2^*  > 


(17c) 


4 


where 


-r-TX  ^  -r^ — j-  \9Q,4.h  -f  —h  sin<^ 

dujdil  [  7o 


fAri/7 

Iq=  dT)’ exp[iTp{xo,0,9)], 

J-^v/2 

fAv/2 

=  /  dr]' sin T}' exp  [ij^{xo,§^,Tj')], 

J—Av/2 


=/■ 


dr}' cos{2t)') exp  [irf>{xo,3,T}')], 


9o,4>  =  fizo  sin  <f>  -  ftys  cos  <f). 


(20a) 


(20fc) 


(20c) 


9o,0  =  {0z(j  cos 4>  +  fiyo  sin  <i>)  cos  ^  ( /3i  -  j  sin 6,  (21a) 

\  ^To  / 


(216) 


Now,  we  will  integrate  over  t}  by  expanding  the  exponential  of  sin  rj  and  cos  rj  in  terms 


of  Bessel  functions 


exp  (td,  sin  2ti]  =  J„(d,)  exp  [i2m77], 


exp  [xd*  cos  17]  =  jr„(d* )  exp  [tn(7r/2  +  t;)]  =  ^  t”  J„(d* )  exp  [inr;].  (23) 

n=— 00  n=  — 00 

Substitute  Eqs.  (22)  and  (23)  into  Eqs.  (20a)-(20c),  we  obtain 


Jo  =  2e'^»  E  i'Mdz)pe,m, 


(24a) 


=  -e‘*"  -,(<«.)  + 


(246) 


(24c) 


Pi,m  —  „• 


2i{2m  +  £  +  do) 


^g»((2»*»+/+<io)Ai7/2]  _  g  — «[(2m+/4-rfo)no/2]^ 


(24d) 
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III.  Fundamental  and  Harmonics 


The  emission  on  axis  is  peaked  at  the  fundamental  and  harmonic  frequencies,  which 
can  be  obtained  from 

Vl,m  =  TtNo - , 

X 

where  x  =  (2m  +  ^  +  do)(7rJVo).  Since  pi^rn  is  peaked  at  x  =  0-  The  frequencies  of  the  on 
axis  radiation  associated  with  the  peak  intensity  are 

l+«V2  +  ^W.+^^.)  ’ 

where  h  ~  —2m  —  £  is  the  harmonic  number. 

The  expressions  for  /o,  /*  and  written  in  terms  of  the  harmonic  number,  become 


Jo  =  2e‘*"  iVl  )A+,„(<i.) 


A=1  m 


~  e  ^  ^  i  Ph  ^  1)  •^m(dx)[Jfc+2m-l(d*)  +  Jfc+2m+l(de)] 

h=l  m 

oo 

Ig  =  — c‘^°  y~^(  — l)”*«Im(d,)[J^42m-2(<^g)  +  •^fc+2m+2(‘^*)), 


(26a) 


(266) 


(26c) 


Pfc  =  ’T-ZVo - , 

Xh 


(26rf) 


where  Xfc  =  (do  -  h)irNo. 

The  analytical  expression  for  the  radiated  energy  per  unit  solid  angle  per  unit  fre¬ 
quency  per  electron,  given  by  the  sum  of  the  expressions  (18)  and  (19),  with  definitions 
given  by  (17a-c),  (21a-b)  and  (26a  d),  are  valid  for  a  large  range  of  values  of  laser  ampli¬ 
tudes.  For  a  <<  1,  only  fundamental  radiation  will  be  observed.  Intensity  of  harmonic 
radiation  becomes  important  for  a  >  1. 

The  spectrum  width  of  a  single  radiating  electron  is 

No 

obtained  by  equating  Xh  =  tt-  Since  the  laser  pulse  typicedly  has  a  large  number  of  periods, 
the  spectral  width  of  the  radiation,  in  principle,  could  be  very  narrow.  However,  the 
spectral  width  will  be  determined  by  the  energy  spread  and  the  emittance  of  the  electron 
beams. 
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rV.  Numerical  Results  and  Comments 


In  this  section  we  present  numerical  results  for  the  energy  radiated  per  unit  solid  angle 
per  unit  frequency  per  electron.  Two  examples  are  given:  1)  no  transverse  beam  velocity 
and  2)  transverse  velocity  /3io  =  0.0075.  In  both  cases,  the  electron  beam  has  70  =  80. 
The  normalized  laser  amplitude  is  a  =  3.2  x  10“^  and  the  number  of  periods  in  the  laser 
pulse  is  taken  to  be  20.  The  small  number  of  periods  is  not  typical  for  a  laser  pulse,  but 
it  illustrates  the  principles,  while  avoiding  difficulties  in  displaying  data  with  very  narrow 
line  widths. 

For  the  electron  beam  without  initial  transverse  velocity,  the  backscattered  radiation 
is  peaked  on  axis.  Figure  2  is  a  plot  of  the  energy  radiated  as  a  function  of  frequency 
and  angle  $  (evaluated  in  the  ^  =  0  plane),  where  wj  is  the  frequency  of  the  fundamental 
based  on  Eq.  (25).  The  bandwidth  should  be  5%,  and  this  is  confirmed  by  Fig.  3,  which 
is  obtained  on  axis.  The  solid  curve  in  Fig.  4  is  a  plot  of  the  peak  energy  radiated  for  each 
given  angle  6.  As  expected  the  radiation  is  confined  within  an  angle  less  than  I/70.  The 
dashed  curve  in  Fig.  4  is  the  corresponding  frequency. 

Realistic  electron  beams  have  finite  emittance.  If  the  electron  beam  has  an  initial 
transverse  velocity,  the  radiation  pattern  is  distorted.  Figure  5  is  a  plot  of  the  energy 
radiated  as  a  function  of  frequency  and  angle  6  (evaluated  in  the  ^  =  0  plane),  for  = 
0.0075.  A  plot  of  the  energy  radiated  on  axis  as  a  function  of  the  normalized  frequency, 
shown  in  Fig.  6,  shows  that  the  radiation  field  on  axis  is  reduced  by  about  25%  as  predicted 
by  the  analytical  expression.  The  solid  curve  in  Fig.  7  is  the  plot  of  peak  energy  radiated 
for  each  given  angle  0  and  the  dashed  curve  in  Fig.  7  is  the  corresponding  frequency.  The 
peak  frequency  is  at  angle  9  ~  while  the  peak  of  the  intensity  is  at  a  smaller  angle. 

We  have  derived  an  analytic  expression  for  the  intensity  distribution  of  Thomson 
scattered  radiation  for  the  case  of  a  linearly  polarized  laser  pulse  incident  on  a  coun- 
terpropagating  electron  beam.  We  have  calculated  the  effects  of  small  initial  transverse 
momentum,  including  the  distortion  of  the  intensity  distribution,  the  reduction  in  on-axis 
intensity,  and  the  increase  in  bandwidth.  We  are  currently  extending  this  work  to  con¬ 
sider  the  effects  of  emittance  (for  a  distribution  of  electrons)  and  laser  pulse  shape  on  the 
scattered  radiation. 
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Fig.  1.  Schematic  of  the  Thomson  backscatter  configi 
sects  an  incoming  laser  pulse. 
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Fig.  2.  Plot  of  normalized  energy  radiated  per  unit  solid  angle  per  unit  frequency  per  electron 
as  a  function  of  normalized  frequency  and  angle  6  (evaluated  in  the  <^  =  0  plane)  for 
an  electron  with  no  initial  transverse  velocity. 


II 


d^I/dwc5Q  (normalized) 
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Fig.  3.  Plot  of  normalized  energy  radiated  on  axis  per  unit  solid  angle  per  unit  frequency  per 
electron  as  a  function  of  normalized  frequency. 
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Q  (normalized; 
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angle  6  (rad) 

Fig.  4.  Plot  of  the  peak  energy  radiated  (solid  curve)  and  the  corresponding  frequency  (dashed 
curve)  for  each  given  angle  9  (evaluated  in  the  <f>  =  0  plane). 
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Fig.  5.  Plot  of  normalized  energy  radiated  per  unit  solid  angle  per  unit  frequency  per  electron 
as  a  function  of  normalized  frequency  and  angle  9  (evaluated  in  the  ^  =  0  plane)  for 
an  electron  with  /9c0  =  0.0075c. 
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Fig.  6.  Plot  of  normalized  energy  radiated  on  axis  per  unit  solid  angle  per  unit  frequency  per 
electron  as  a  function  of  normalized  frequency. 
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Fig.  7.  Plot  of  the  peak  energy  radiated  (solid  curve)  and  the  corresponding  frequency  (dashed 
curve)  for  each  given  angle  9  (evaluated  in  the  ^  =  0  plane). 
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